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Abstract 

With the development of high-strength NdFeB magnetic material, moving magnet linear compressors have been gradually introduced 
in the fields of refrigeration and cryogenic engineering, especially in Stirling and pulse tube cryocoolers. This paper presents simulation 
and experimental investigations on the static and dynamic characteristics of a moving magnet linear motor and a moving magnet linear 
compressor. Both equivalent magnetic circuits and finite element approaches have been used to model the moving magnet linear motor. 
Subsequently, the force and equilibrium characteristics of the linear motor have been predicted and verified by detailed static experimen¬ 
tal analyses. In combination with a harmonic analysis, experimental investigations were conducted on a prototype of a moving magnet 
linear compressor. A voltage-stroke relationship, the effect of charging pressure on the performance and dynamic frequency response 
characteristics are investigated. Finally, the method to identify optimal points of the linear compressor has been described, which is indis¬ 
pensable to the design and operation of moving magnet linear compressors. 

© 2007 Published by Elsevier Ltd. 
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1. Introduction 

Linear free piston compressors are employed in a wide 
range of reciprocating electro-mechanical actuation sys¬ 
tems, from domestic refrigerators, Stirling cycle and pulse 
tube cryocoolers, to fuel pumps and artificial heart devices. 
They are usually required to be efficient and reliable, typi¬ 
cally offering 20,000 h of maintenance free operation. 
Moreover, they must also meet demands for functionality 
and cost, whilst fulfilling minimum volume and loss criteria 
[1]. Linear motors which are the heart of the linear com¬ 
pressor can be of three different types, i.e. moving coil, 
moving magnet and moving iron. The moving iron type 
is scarcely used due to the heavy moving parts and small 
force. The moving magnet and moving coil motors have 
their own advantages and disadvantages. Flying leads in 
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the moving coil motor is a potential source of unreliability 
and limits the achievable stroke. On the contrast, due to the 
windings bonded directly to a yoke in the moving magnet 
motor, the good thermal dissipation characteristics, low 
material outgassing rate and high reliability can be 
obtained. Moreover, the required volume of the permanent 
magnet (PM) in the moving magnet configuration is much 
smaller than that needed by the moving coil configuration 
of same power output and efficiency. All above character¬ 
istics make the moving magnet configuration more suitable 
to be used in the linear compressor [2]. Recently, the 
increasing number of moving magnet linear compressors 
has been used in some important space missions [3,4]. 

Redlich et al. [1] introduced a type of linear motor, 
whose shuttle consisted of a ferromagnetic core and a 
homopolar, radially magnetized, bonded PM ring, a stator 
comprised a pair of coils, wound in series opposition and 
an outer ferromagnetic sleeve (see Fig. 1). This topology 
has been widely used in gas compressors. From the 
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Nomenclature 



B 

magnetic density, T 

</> 

magnetic flux, Wb 

D 

damp, Ns/m 

•A 

flux linkage, Wb Turn 

e 

back voltage, V 

CO 

angular frequency, rad/s 

F 

force, N 



H 

magnetic potential, A/m 

Superscripts 

i 

current, A 

a 

amplitude 

I 

effective value of current, A 

c 

coercive 

j 

imaginary part of complex 

Cu 

copper 

K 

stiffness, N/m 

e 

equivalent 

i 

length, mm 

g 

gas 

L 

inductance, H 

in 

input 

M 

mass, kg 

iron 

yoke iron 

N 

number of windings 

m 

mechanical 

P 

power, W 

max 

maximum 

P 

pressure, Pa 

mean 

mean value 

r 

radius, mm 

mover 

piston and actuator of motor 

R 

electrical resistance, Q 

n 

natural frequency 

S 

area, m 2 

0 

operation 

t 

time, s 

P 

permanent magnet (PM) 

T 

period, s 

pi 

left PM 

V 

voltage, V 

pr 

right PM 

V 

volume, m 3 

r 

remnant magnet 

w 

co-energy, J 

s 

spring 

X 

displacement from balance point, mm 

t 

piston 



V 

virtual 

Greek symbols 

w 

wall 

n 

efficiency, % 

0 

vacuum 

e 

phase angle, rad 



fi 

relative magnetic permeability 




viewpoint of electro-mechanical energy conversion, the 
electrical energy in the exciting circuit is converted into 
mechanical energy by means of induced electromotive force 
(EMF) produced by the periodical variation of magnetic 
flux in an air gap. So the distribution and organization of 
the magnetic field in the air gap have great effects on the 
motor performance. Until now, the equivalent magnetic 
circuit analysis and finite element analysis (FEA) are two 
of the most effective methods in the magnetic field analysis. 


h 



Fig. 1. Definitions of the geometrical dimensions of the moving magnet 
linear motor. 


For the equivalent magnetic circuit analysis, Fongthworth 
[5] gave the simplified expression of force and back EMF, 
and the engineering design and construction method of a 
dual PM linear motor in a Stirling cryocooler was intro¬ 
duced. Byckling and Perkio [6] developed the equations 
of motion, force and the equivalent circuit of a moving 
magnet transducer for a loudspeaker. It was pointed out 
that the moving magnet transducer had low pass and 
high-frequency cut-off characteristics. Yang et al. [7] gave 
the analytical solution of the electromagnetic force of a lin¬ 
ear oscillating motor and simulated the field, force and 
dynamics respectively. Poornima and Hsu [8] gave the 
expressions of force and the electrical power converted to 
produce motion in a moving magnet loudspeaker system. 
Ebihara and Watada [9] discussed the static thrust and fre¬ 
quency characteristics of a linear oscillatory actuator used 
in an artificial heart. The maximum static thrust was mea¬ 
sured at approximately 72 N and the starting thrust at 
approximately 20 N. Clark et al. [10,11] presented an opti¬ 
mum design method of the linear motor based on the FEA. 
The Coil inductance and peak force/amp of their prototype 
actuator were 1.65 mH and 8.0 N/A, respectively. These 
investigations show that the FEA is an effective method 
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to get the distribution parameters of the magnetic field in 
the linear motor. Although the variation of force along 
the stroke was actually observed in the experiments 
[9,10,12], the electromagnetic force along the whole stroke 
was still regarded as a constant in majority of above-men¬ 
tioned papers. For the practical application of the moving 
magnet linear compressors, this characteristic is very 
important owing to the increasing gas stiffness in the pro¬ 
cess of compression. 

Another two important aspects of the moving magnet 
linear compressor are the thermodynamic process of work¬ 
ing gas and the dynamic process of the piston, both of 
which are similar to the moving coil linear compressor. 
Many previous papers [13-17] reported the dynamic char¬ 
acteristics of the moving coil compressors, in which the 
effect of the compressed gas on the piston had been 
expressed by an equivalent damping coefficient and a gas 
spring stiffness. The simulating results of the actuator are 
in good agreement with experimental data. 

Although Sunpower Inc. produced a series of high per¬ 
formance Stirling and pulse tube cryocoolers based on the 
moving magnet linear compressors, almost all published 
papers [18-20] introduced the performance as a whole 
rather than the detail analysis of the parameters. Philips 
Laboratories [21] provided a simplified design method of 
the moving magnet linear compressor for the Stirling 
cooler. However, the variation of the electromagnetic force 
and inductance of the coil were ignored. 

In this paper, the force of a linear motor, as the function 
of geometrical parameters, input current and magnetic 
potential of the PM is firstly deduced from the magnetic cir¬ 
cuit analysis. Then, the FEA was used to obtain the distri¬ 
bution of the magnetic field in the linear motor. Based on 
the experimental measurements, the static performance of 
the linear motor, including force and equilibrium character¬ 
istics have been studied. Furthermore, in order to 
investigate the performance of a linear compressor compre¬ 
hensively, a mathematical model has been developed, in 
which the static force model, thermodynamic and electrical 
circuit model have been involved. Based on harmonic 
assumptions, the two dimensional control equations of lin¬ 
ear compressors have been solved and consequently the 
dynamic characteristics of the linear compressor have been 
revealed by a series of equations. As the result, the corre¬ 
sponding equivalent electrical circuit of a linear compressor 
is presented, by which the compressor can be regarded as a 
simple input and output module. Finally, a prototype of the 
moving magnet linear compressor was constructed and 
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Fig. 2. Magnetic equivalent circuit of the moving magnet linear motor. 

tested. Theoretical solutions have been verified by the 
experimental data and many useful results are presented. 


2. Moving magnet linear motors 

2.7. Magnetic circuit model of linear motors 


The electromagnetic force of linear motor depends on 
the interaction field between the exciting coil and the 
PM. Referring to the theory of electromagnetic analogy, 
the magnetic circuit of the moving magnet linear motor 
can be simplified as shown in Fig. 2, where the linear 
B-H characteristics of the PM and the lumped approxima¬ 
tion of the distributed flux are adopted [6]. The quantity <A 
is the flux in the air gap which goes through the permanent 
magnet and that of </>2 is the flux that bypass it. The quan¬ 
tity Hq is the coercive force of the PM, and N and / are the 
number of turns and current of the driving coil. Due to lin¬ 
ear assumptions of the magnetic circuit, the magnetic co¬ 
energy of this system in Eq. (1) can be divided into three 
parts which correspond to three sources of the magnetic 
potential in Fig. 2 respectively 


W(iy pi, Apr,*?*) — / dfF(z, z'vpl — Apr — 0,x) 

J Path 1 

+ / dfT(Apr, A P i = 0, i = fx) 

J Path2 

+ / d W (Apl f = f ,Apr = 7vpr A) 7 
J Path3 


(i) 


where Pathl, Path2 and Path3 indicate three different inte¬ 
grating paths to get the magnetic co-energy (fF(/ Vp i, 
Apr? 7))* The variables of Api and Apr are respectively de¬ 
fined as the virtual currents, which can bring up the same 
magnetic potential of driving coil as that of the PM [22]. 
The quantity of / Vpr and / Vp i are their corresponding values 
of NdFeB. Based on Ohm’s and Kirchhoffs law of the 
magnetic circuit, the energy stored in the system is x 
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where R\, R 2 , R 3 and R 4 are the reluctances of the air gap 
corresponding to the magnetic flux </> 1? 0 2 , <£3 and 0 4 , 
respectively; and i ^ 2 are the reluctances of the PMs, 
respectively; Ri ron is the reluctance of the iron core. 

In the same way, Hp a th 2 and Hp a th 3 can be obtained by 
integration of their own paths, as follows: 

iPi+R 2 + R non )(H 0 l p ) 2 /2 - R 2 NIH 0 l p 

RliP\ + Rnon ) + (^1 + Rpi) (p | + Rl +^iron) 

and 


IP Path3 


^R 2 NIH 0 l p + 


(go/p) 

M 



(q l -q 2 )+q l R m 


q 2 R 4 , 



where the intermediate variables are 


Pi = ^ 4(^3 + Rpl) / {R 3 + + ^ 2)5 

#1 = R 2 R 4 + ^ 4(^1 + 

#2 = ^2^iron + (f^2 + ^iron)(^l + 

and 


one can get / =-^^-H 0 l p (H 0 l p fR). When NI = 

—H 0 l p = R(j)/ 2, the maximum force can be obtained. How¬ 
ever it should be pointed out that it is impossible to raise 
NI infinitely by increasing the input current in reference 
to the limit varying range of H o (0,H c ). So, the maximum 
permitted current (/ max ) is confined by H c l p = M max . In 
other words, the maximum permitted peak current is 
/max = U c l p /N so that the maximum force is / max = 
2nrn Q H 0 H c l 2 v /l\. This expression is the theoretical value 
of the force. In order to make full use of the PM, the oper¬ 
ation point of the PM should vary closely to the point with 
the maximum magnetism energy produce. As to the 
NdFeB, the coercive force of this point is about 0 .5H C . 
Therefore, in the design of linear motors, both magnetic 
circuit organization and property of PM should be consid¬ 
ered comprehensively. 

According to the Faraday’s law of the electromagnetic 
induction, the induced back EMF applied on the exciting 
coil can be expressed as 

0i/A di /0t/A dx 
0 /) dt \ 0x ) dt 1 




M — [R 2 R 4 +^ 2(^3 T Rp2) \ [(^2 +^iron)(^l + Rp \) + f^2^iron] 
+ ^2(^3 + R p 2 )[ R 2 Ra + ^4(^1 +^l)]- 

Then, the electromagnetic force equals the partially deriva¬ 
tive of the IF in Eq. (1) with respect to the displacement (x), 
that is 


where the W is the coil flux linkage. 

Referring to Fig. 2, the induced back EMF on the excit¬ 
ing coil can be obtained by the magnetic circuit analysis. 
That is 


e = N 2 


nrji^a 


di 

dt 


~NH 0 l p 



dx 

dt' 
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dW 



z Vp 14 Vpr 4 



Substituting the result of Eq. (1) into Eq. (5), the electro¬ 
magnetic force is 
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InrjiQ 2nrf.i 0 x . 2 2nrfi 0 2x 
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H 0 IpNI - 


h a {Holp) + h 


h~l 


m\ 



By means of dimensional analysis, it is found that the 
quantity of the inductive voltage (the first term of Eq. 
( 8 )) is far less than that of the speed voltage (the second 
term of Eq. ( 8 )). So, the majority of the electrical energy 
is actually converted into the mechanical energy by the 
virtue of the speed voltage. 

2.2. Static performance of linear motors 


which indicates that the driving current, geometrical 
parameters and characteristics of the PM have close rela¬ 
tionship with the force. The first term actually represents 
interacting force between the PM and the coil. The last 
two terms are the functions of the displacement and be¬ 
come zero in the central position. By the dimensional anal¬ 
ysis, it was found that the absolute value of the first part is 
larger than that of the last two parts. 

In practice, the operating point of the PM actually var¬ 
ies along the magnetization line. H 0 in Eq. ( 6 ) is the mag¬ 
netic potential of the operating point instead of the 
coercive force (H c ). As for the NdFeB, the magnetization 
line in the second quadrant can be approximately expressed 
as a linear function, B = —HB r /H c + B y . The operation 
point along the magnetization line is determined by the flux 
across the PM. Assuming the summary of the magnetic 
potential provided by the PM is H 0 l p , according to the con¬ 
tinuity of the magnetic potential (NI = H 0 l p + fR), and 
substituting this expression into the first term of Eq. ( 6 ), 


The static force-displacement characteristic, i.e. the rela¬ 
tionship between the electromagnetic force exerted on the 
moving PM and its axial displacement under different 
direct currents applied on the exciting coil, is one of key 
indexes to the linear motor. The magnetic field distribution 
and the magnetic force without exciting current also have 
great influence on linear motor applications. 

The parameters of the geometrical dimension, electrical 
configuration and material property are given in Table 1. 
The diameter of the investigated linear motor stator (see 
Fig. 3) is 120 mm, the length is 65 mm and the mass is 
2.21 kg. In order to reduce the Joule heat (I~R), the 
large-diameter copper wire was wound tightly on the yoke. 
The total electric resistance of the coil is 0.41 Q. The mag¬ 
netic circuit is formed by ferrite, which can effectively 
reduce the core loss. The diameter of the corresponding 
aluminum armature, attached by the lamellar NdFeB-35 
PM, is 72 mm and its mass is 0.41 kg. The static character¬ 
istics of this motor were investigated on a modified C620H 
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Table 1 

Parameters of the moving magnet linear compressor 
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Linear motor 

Magnet material 

Relative permeability 

1.0998 

Coil and yoke 

Outer diameter (mm) 

94 

Coercivity (A/m) 

-8.9 x 10 5 

Inner diameter (mm) 

76 

Retentivity (Tesla) 

1.23 

Coil wire size (mm) 

1.5 

Thickness (mm) 

4 

Number of turns 

150 

Magnetization direction 

Actuator 

Radial 

Coil resistance (H) 

0.41 

Length (mm) 

45 

Max. current (A) 

15 

Diameter (mm) 

73 

Max. electrical input (W) 

250 

Mass (g) 

410 

Distance between internal and external yoke (mm) 

6 

Integral linear compressor 

Frequency (Hz) 

45-65 

Total stiffness of spring (N/m) 

2316 

Max. displacement (mm) 

10 

Total mass of spring (g) 

622 

Diameter of cylinder (mm) 

45 

Total mass of mover (g) 

560 

Volume of gas bank (cm ) 

141 

Clearance seal (pm) 

10-20 

Working gas 

Helium 





Fig. 3. Stator of the linear motor. 


machine tool (see Fig. 4), on which the force and displace¬ 
ment could be conveniently regulated and kept coaxial. An 
ergometer (TK-10) and a slide calliper were used to mea¬ 
sure the force and displacement respectively. In addition, 
magnetic field analysis software (Ansoft Maxwell) was used 
to get the optimum distribution of the magnetic field. By 
adopting different profiles of the yoke, even distribution 
of flux in the air gap had been obtained, at the same time 
the maximum value of flux was kept less than the magnetic 
saturation point of the ferrite. 

Fig. 5 presents the FEA results of the magnetic filed dis¬ 
tribution in the linear motor. It was found that the mag¬ 
netic field in the air gap is the superposed fields of the 
PM and the exciting coil. From the density difference of 
flux at two ends of the moving PM (Fig. 5a), it is evident 


that the force directs upward and the subsequent motion 
is to offset this difference. The flux lines are not parallel 
with each other in the air gap, especially at the both ends 
of the PM. It means that there exist axial and radial com¬ 
ponents of the force. The axial component is the useful part 
and the radial component counteracts circumferentially. 
Due to the proper section configuration of the yoke, the 
eddy current is actually not obvious. The maximum flux 
density can be kept less than 2.0 T and the magnetic satu¬ 
ration can be avoided. However, there exists weak leakage 
of flux at the end of the yoke because of the limited 
dimensions. 

According to the equipotential line of the magnetic 
potential (see Fig. 5b), it was found that the operating 
potential of the PM is actually far less than the coercive 
force of NdFeB-35. In the practical operation, the flux den¬ 
sity across the PM rarely declines to zero. Based on Fig. 5b, 
the practical operating point of the PM can be approxi¬ 
mately obtained. 

In a moving coil linear motor, the magnetic force in the 
air gap depends on the existence of exciting current. On the 
contrast, the static force equilibrium of the PM in the mov¬ 
ing magnet linear motor is more complex. Without exciting 
current, the magnetic force is still applied on the PM. In 
order to investigate this special phenomenon, both an 
experimental study and a FEA have been conducted. 

In the experimental study, the magnetic force was mea¬ 
sured when the stroke of the PM range from 0 mm to 
9.0 mm. As shown in Fig. 6, when the PM is in the central 



Fig. 4. Schematic of the static characteristic testing bench. 
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Flux Lines [Wb] 



2.1635e-004 
2.0135e-004 
1.8634e-004 
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1.4133e-004 
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2.129Se-005 
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-8.7132e-006 
-2.3717e-00S 
-3.8722e-005 
-5.3726e-00S 
-6.8730e-00S 








H[A/m] 

1.1212e+006 
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8.2613e+005 
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5.9009e+004 
0.OOOOe+OOO 


Fig. 5. The magnetic field distribution at the exciting current of 9 A. (a) Magnetic flux, (b) Magnetic potential. 



Displacement from the Central Point(mm) 


Fig. 6. Magnetic force without exciting current with displacement of the 
shuttle. 


position, the composite magnetic force applied on the PM 
is equal to zero. However, it was observed that once the 
disturbance is applied, the magnetic force will increase rap¬ 
idly with the increase of displacement. In the experiments, 
the increasing rate is about 4 N/mm. As a result, the equi¬ 
librium will be broken easily. Hence, this state can be 
regarded as the unsteady equilibrium. 

This magnetic force characteristic without the exciting 
current should be considered in advance in the design of 
the linear motor. Once the PM leaves the central point 
due to the inevitable outside distribution, enough spring 
stiffness should be provided to avoid large offset displace¬ 
ment. Under this off-center condition, more input current 
should be provided to ensure the oscillation around the 
central point. It also means the enough input current 
should be provided to ensure that the first term in Eq. (6) 
is larger than the last two terms. 

When the displacement increases to 11 mm continu¬ 
ously, the magnetic force correspondingly reaches the max¬ 
imum value of 45 N according to the FEA results. Then, 
this force begins to reduce gradually and finally becomes 
zero when the displacement is about 26 mm. Similar to that 


of the central point, the composite force at this point is also 
equal to zero. However, the essential difference between 
these two points is that the magnetic force is negative 
(see Fig. 6) when the displacement is larger than 26 mm. 
It means that once the PM leaves this point, the magnetic 
force will draw it back. So the state of this point can be 
considered as a steady equilibrium. The position of this 
point in the experimental study was 24.4 mm, which is close 
to the result of the FEA. 

The force characteristics of these two equilibrium points 
imply that, due to the inevitable disturbance, the piston of 
the linear compressor actually deflects its central point, 
which results in a reduction of the compressing stroke. 
Therefore, this factor should be considered in the design 
of mechanical springs to avoid the excess deflection. This 
force has close relation to the geometrical dimension of 
the yoke and PM. According to Eq. (6), if the dimension 
of the air gap and the diameter of the yoke are constant, 
the force should be the function of 2 /(h — l) — 1 /a. When 
2a = h — /, referencing to Fig. 1, this force will tend to zero. 
Hence, this force can be reduced significantly by increasing 
the length of the yoke. 

Using three different methods, i.e. experimental mea¬ 
surement, theoretical calculation from Eq. (6) and FEA 
respectively, the relation between the axial displacement 
and the force at the exciting current of 11 A is illustrated 
in Fig. 7. It is shown that the electromagnetic force declines 
with the increase of the displacement. For example, the 
force declines about 23.8% with the displacement varying 
from 0 mm to 8 mm. It can also be found that both the 
FEA and linear magnetic circuit analysis are the effective 
approaches to study the force-displacement characteristics 
of the linear motor. Due to the existence of the flux leakage 
at the end of the stroke, which had not involved in the mag¬ 
netic circuit model, there is about 20% discrepancy at full 
stroke. 

It was found from Fig. 8 that the electromagnetic force 
at the central position increases when the exciting current 
increases from 6 A to 11 A. When the PM is at the central 
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Displacement from the Central Point(mm) 


Fig. 7. Variation of the electromagnetic force with the displacement. 


and the equivalent damping of the compressing gas (Z> e ); K 
is the sum of the flexure spring stiffness (K m ) and the equiv¬ 
alent stiffness of the gas spring (K g ). The gas spring effect 
on the piston can be expressed as the equivalent gas stiff¬ 
ness as follows: 

K g =—————————————. (ip 

X 

The equivalent damping coefficient (Z> e ) of the compressed 
gas, reflecting the hysteresis loss caused by irreversible 
compression, as well as the viscous damping coefficient 
(Z) g ) can be obtained respectively in references [23,24]. 

3.2. Harmonic analysis 



Fig. 8. Variation of the magnetic force with the input current. 


position and the current increases from 6 A to 9 A and 9 A 
to 11 A, the corresponding increment rate of the force are 
44.06% and 23.08% respectively. It can also be observed 
that the input current is approximated proportional to 
the force when the actuator is located at the central point. 
The comparison between the results from the FEA and Eq. 
(6) with the experimental data shows the validity of the the¬ 
oretical analysis to predict the force characteristics of the 
linear motor. 

3. Theoretical analysis of linear compressors 

3.1. Mathematical model of linear compressors 

The voltage equation for the driving circuit of a linear 
motor can be expressed as 

v(t) = iR + e. (9) 

The mechanical force balance of the moving PM screwed 
on the piston results in the following equation 

d 2 x dx . . 

M df +D d~ t + Kx =f(x,i), ( 10 ) 


In order to obtain explicit solutions of Eqs. (9) and (10), 
the harmonic wave assumptions of input voltage, electrical 
current and displacement of the piston were used as 

x(t) =x a e j ^, (12) 

i(t) =/ a e j( ^ +0l) , (13) 


and 

v (t) = v^ cot+62 \ (14) 

Although there exist higher-order harmonic waves, the 
deformation of sine wave caused by them is actually small, 
and penalty of correction can be ignored [3,17]. As a result, 
the relation between the piston displacement and the cur¬ 
rent is 


*a _ _ (h _ 

h f{A\ - M ® 2 ) 2 + (£><u ) 2 



and the phase shift between the current and the displace¬ 
ment is 


0 \ — arctan 


Dm 


A x - Men 2 




From Eqs. (15) and (16), it was found that when A\ equals 
to Mco or in other words, when the resonant state of oscil¬ 
lation system is reached, the largest amplitude of the dis¬ 
placement can be obtained with a certain current input. 
In the same way, the relation between the piston displace¬ 
ment and the voltage is 


X, 


V, 


f(A x -Mco 2 )- 


fDoj 2 

a 2 


i 2 


+ 


coiR^+A.+f^-Mco 2 ) 


1 2 


(17) 


and the phase shift between the applying voltage and dis¬ 
placement is 


0 2 = arctan 



(Ry 2 +A A +f i {A l -Mco 2 )) 

£04,-Mco2)-fDco2 



where M is the total mass of the moving PM and piston; D As a result, the phase shift between applied voltage and in¬ 
is the sum of the viscous damping of the clearance seal (Z) g ) put current is 
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0 = 0] — Oj = — arctan 


(A 3 D 2 co 3 -\-coA 2 A 4 (Ai —Mco 2 ) +cdA 3 (Ai — Mco 2 ) 2 ) 


R(A i -Ma > 2 ) 2 + a> 1 D(RD+A 2 A i ) 


(19) 


Substituting the expression of the electromagnetic force 
(Eq. (6)) into the following equation and integrating, the 
mechanical power of the linear motor is 

1 ^ 

Pm = T I /(*>*) &C = c0 (° 7ir H 0 l v NxJz sin Q x 
1 Jo n 

0.5 Dil 

_ Cl 

f A | -Mo) 2 \ 2 , 

V ^ 



The electrical input power of the linear motor is 

P -Al, ' d ' = 2 

and the copper loss of the motor is 


>r 1 
vdi = -z a r a cos 0, 


( 21 ) 


Pcu = J fRdt = ^i 2 d R. 


( 22 ) 


Finally, the motor efficiency defined as the ratio of the 
mechanical output from the motor to the electrical input 
power is 
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p 


m 


Pi 


m 


R 


' d\ I ( a { -MoA 
a 2 ) " r ^ A 2 (D J 


(23) 


ED 


It can be seen that when the input power (p m ) is constant, 
the motor efficiency depends on two groups of design 
parameters. The first group is related to the operating con¬ 
ditions. The motor efficiency can be maximized when reso¬ 
nance is realized. The second group is determined by the 
motor size and material properties. The motor efficiency 
can be improved greatly when (1) the air gap become small, 
(2) both the magnet and coil sizes are large enough and (3) 
the coercive force of the PM becomes large. 


3.3. Equivalent circuit of linear compressor 

The compressor can be regarded as an inductance com¬ 
ponent which is important to the design of the drive and 
control circuit [24]. By combining Eq. (9) with Eq. (15), 
we have 

t = (RRe)j(LL g ), (24) 

l 

where 

DA 2 or 

n _ _ l _ 

e ~~ (Ai — Mco 2 ) 2 + (Deo ) 2 ’ 
co(Ai — Mco 2 )A 2 
6 ~~ (A, - Mco 2 ) 2 + {Deo) 2 ' 


Here R is the equivalent resistance of the copper loss and 
iron loss. R e is the equivalent resistance of the mechanical 
damping. L is the inductance of the coil and L e is the equiv¬ 
alent inductance of the mechanical springs. 

As a consequence, the moving magnet linear compressor 
can be expressed as an equivalent circuit, as shown in 
Fig. 9. According to the electric circuit theory, the induc¬ 
tive components absorb and discharge energy periodically, 
and the energy dissipation just happens on the resistance. 
Similarly, the mechanical power of the linear compressor 
can be regarded as the energy dissipation on the R e , so it 
can be expressed as 0.5z^ e - The motor efficiency can finally 
expressed by 

n =-. 

1 r+r q 

Substituting the expression of R e and R into Eq. (25), the 
motor efficiency can be obtained which is identical to Eq. 
(23). 

4. Experimental study of linear compressors 

In order to investigate the dynamic performance of the 
linear compressor, a measurement system was built (see 
Fig. 10). The above-mentioned linear motor was used to 
provide driving force for compression. The moving magnet 
shuttle was supported on stainless steel suspension flexure 
springs, whose stiffness is shown in Table 1. A hollow pis¬ 
ton was adopted to reduce the mass of the reciprocating 
assembly. In order to investigate the effect of piston mass 
on the natural frequency, two pistons with the identical 
dimensions were manufactured by aluminum (165 g) and 
stainless steel (483 g) respectively. The clearance seal 
between the piston and cylinder is 10-20 pm in width and 
80 mm in length. The static leakage ratio to the swept vol¬ 
ume is consequently less than 5%. A pure sine-wave power 
(TPS-500B) supplies single-phase alternating current to the 
primary coil of transformer (TDGC2M-5, 5 kW), the 
motor exciting coil has been connected with the secondary 
coil serially. As a result, large exciting current (maximum 
15 A) can be provided to produce enough magnetic flux. 
The dynamic pressure signal was measured by the piezo¬ 
electric transducers (KISTEL 601 A) and recorded by a 
digital oscilloscope. The displacement of the piston was 
measured by the strain signal of the flexure spring. Table 2 
is the list of adjustable working conditions, by combination 
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Fig. 9. Equivalent circuit of the linear compressor. 
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Fig. 10. Experimental apparatus of the moving magnet linear compressor. 


Table 2 

Working conditions of the linear compressor 

Charging pressure (MPa) 0.8 1.0 1.2 1.4 1.6 1.8 

Input voltage (V) 5 10 15 20 24 

Working frequency (Hz) 45 50 55 60 65 


of which different experimental conditions had been 
realized. 

Fig. 11 shows the influence of input voltage on the dis¬ 
placement of the piston under different operating frequen¬ 
cies. It is evident that the relationship between the input 
voltage and displacement is approximately linear and the 
increasing rates vary with the working frequency. At the 
working frequencies of 45 Hz and 65 Hz, the linear increas¬ 
ing rates are 0.47 mm/V and 0.12 mm/V, respectively. 
When the working frequency reaches its natural frequency, 
the increasing rate at this point is larger than that at other 
conditions. 

The predicted strokes from Eq. (17) agree with the 
experimental results, especially at high working frequen¬ 
cies. The stiffness of the gas and damp coefficient play a 
key role on the accuracy of the prediction. In practice, 
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Fig. 11. Comparison of predicted voltage-stroke characteristics for 
different operating frequencies, with measured characteristic, at a charge 
pressure of 1.2 MPa. 


the damp mainly involves the mechanical friction damping 
and the thermodynamic dissipation of the compression gas. 
Among them, the mechanical friction, depending on the 
magnitude of the clearance, can be measured in the process 
of assemblage by virtue of the vibrometer. 

Fig. 12 shows the variations of the measured amplitude 
of the displacement, the pressure wave and the correspond¬ 
ing natural frequency under different working frequencies 
with the charging pressure ranging from 0.8 MPa to 
1.8 MPa, at the constant input voltage of 20 Y. The natural 
frequency can be obtained by — \JK/ (M mover + M s /3), 
in which the spring stiffness of gas (K g ) is from Eq. (11). 
It is clear that with the increase of charging pressure, nat¬ 
ural frequency rises. When the working frequency is close 
to its natural frequency, larger stroke can be obtained. 
When the charging pressure is 1.4 MPa, the corresponding 
natural frequency is approximately 46 Hz, and the stroke 
at 65 Hz working frequency is smaller than that at 45 Hz. 

The pressure wave actually does not reach its peak with 
the displacement wave synchronously. If the total mass of 
the gas in the compression chamber is assumed constant, 
the ideal gas law equation can be transformed to A P/ 
A V = Pmean/ F me an ? i n which A V equals to S t Ax. It can be 
found that the amplitude of the pressure wave is related 
to the charging pressure and the variation of the piston 
stroke. 

Finally, the natural frequency is not a simple function of 
charging pressure and is also influenced by the stroke [25]. 
When the working frequencies ranges from 45 Hz to 55 Hz, 
the influence of them on the natural frequency is less than 
10% with the charging pressure ranging from 0.8 MPa to 
1.8 MPa (see Fig. 12). However, when the working fre¬ 
quency is 65 Hz, the corresponding natural frequency 
(40 Hz) is very low due to the weak stiffness of the com¬ 
pressed gas. 

The frequency responses of two different pistons are 
shown in Fig. 13. The mass of the two shuttles are 767 g 
and 1085 g respectively. As for the aluminum piston, at 
the charging pressure of 1.6 MPa and the applied current 
of 10 A, the measured pressure amplitude is 365 kPa, the 
peak-peak stroke is 11.34 mm and the natural frequency 
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0.8 1.0 1.2 1.4 1.6 1.8 



Fig. 12. Measured displacement, pressure wave and natural frequency vs. charging pressures for different operating frequencies at 20 V input voltage. 



Working Frequency(Hz) 


Fig. 13. Comparison of predicted frequency-stroke characteristics for 
different shuttles, with measured characteristics, at a charge pressure of 
1.6 MPa and 10 A input current. 

is 48 Hz. When the exciting current is constant, there exists 
a peak of stroke when the working frequency is close to its 
natural frequency. Comparatively, the natural frequency of 
the steel piston is about 40 Hz at 1.6 MPa charging pres¬ 
sure and 10 A input current. In consequence, the stroke 
decreases with the increase of the working frequency. 
Therefore, the reasonable configurations and materials 
should be adopted to ensure good system performance. 

Neglecting the influence of other parameters in Eqs. (20) 
and (23), both efficiency and mechanical work of the com¬ 
pressor are decided by the term of — Mco 2 )/'A 2 co. It can 
be found that when ^ Mco , this term will decrease with 
the increase of working frequency; when A\ < Moo , the 
value of this term will increase with the increase of working 
frequency. In Fig. 14, when the working frequency 
approaches to its natural frequency, the mechanical work 
is always higher. However, the point of maximum efficiency 
is not the point of maximum work. For instance, at the 



45 50 55 60 65 

Working Frequency(Hz) 

Fig. 14. Measured motor efficiency and mechanical work vs. working 
frequency at the input voltage of 15 V and different charging pressures. 

charging pressure of 1.8 MPa, the natural frequency is 
about 50 Hz. The peak of efficiency 83% happens at the 
working frequency of 45 Hz, contrastively the maximum 
work output of 134 W happens at the resonant point. 
According to Fig. 14, it can be found that the efficiency 
decreases and the output work reaches its peak at the res¬ 
onant point. When the resonant happens, the stroke is cor¬ 
respondingly maximized, and the damping reaches its peak. 
Hence, the efficiency may decrease a little, according to Eq. 
(23). In the experimental study, it can be observed that 
when the resonance happens, the current was always be 
maximized with the same input voltage and the mechanical 
work was also be maximized. In addition, the maximum 
current as well as the stroke will produce large energy dis¬ 
sipation. That is the reason why the maximum efficiency 
cannot obtained at the resonant state. Generally, the work¬ 
ing frequency should be regulated a bit lower than its 
































































N. Chen et al. / Cryogenics 47 (2007) 457-467 


467 


natural frequency. Only at this condition, the relative high 
efficiency and large output work can be obtained 
simultaneously. 

5. Conclusions 

The static force characteristics of the moving magnet 
linear motor have been studied by the magnetic circuit 
analysis and numerical FEA, and verified by experimental 
investigations. The expression of the force against the dis¬ 
placement of the piston was deduced, which can actually 
be divided into three parts, one is the force caused by the 
exciting current; the other two terms whose magnitude 
has close relationship with the displacement of the piston, 
are the force interacting between the PM and the ferrite 
yoke. Based on these expressions, the optimization operat¬ 
ing point of the PM was analyzed. The static equilibrium 
characteristics of the PM, including two steady equilibrium 
and one unsteady equilibrium points, and its effects were 
also revealed. 

In order to investigate the performance of the linear 
compressor, a mathematical model has been developed, 
which includes the model of the linear motor, thermody¬ 
namics of compressing gas and driving electric circuit. 
Based on the harmonic solutions of the mathematical 
model, an equivalent circuit of the linear compressor was 
deduced. A set of expressions concerning the mutual rela¬ 
tions between the voltage, current, stroke, phase shift, 
mechanical work and motor efficiency were obtained. In 
addition, a prototype of the linear compressor was devel¬ 
oped and corresponding experimental analysis were con¬ 
ducted. Both the experimental and theoretical results 
indicate that the motion of the moving piston complies 
with the principle of oscillation dynamics. The stroke is 
approximately a linear function with the input voltage. 
The optimal working frequency should be a bit lower than 
its natural frequency. 
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